J/il) Production Via Photon Fragmentation at Hadron Colliders 
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The transverse momentum (pt) distributions of production and polarization for J/^p measured 
by CDF Collaboration are still challenging our understanding of the heavy quarkonium production 
mechanism even with recent significant theoretical progresses on the next-to-leading-order (NLO) 
QCD calculation. We suggest a new mechanism for J/tp production at hadron collider, pp{p) 
7*(J/V') -|- X with J/tp from a virtual photon 7* fragmentation. Our calculations show it's pt 
distribution at NLO will be larger than that of the conventional J/tp production from the color- 
singlet mechanism at NLO when pt > 26 (35) CeV at the Tevatron (LHC) and reach about 6 
(10) times of the conventional one when pt = 50 (100) GeV at the Tevatron (LHC), in spite of 
a suppression factor (a/as)^ that is associated with the QED and QCD coupling constants. In 
addition, it also has large impact on the pt distribution of J/tp polarization in large pt region. 
Therefore, it is a important mechanism for J/tp production at large pt region especially for the 
LHC. 

PACS numbers: 12.38.Bx, 13.66.Bc, 14.40.Gx 



As good probes for people to understand both the per- 
turbative and non-perturbative aspects of QCD, the pro- 
duction of heavy quarkonium in high energy processes 
has attracted considerable attention in recent years(see|l| 
for reviews). 

Within the framework of an effective field theory non- 
relativistic QCD (NRQCD) 0, the orders of magnitude 
discrepancies between experimental data measured by 
the CDF Collaboration[3j and the color-singlet mecha- 
nism (CSM) theoretical predictions were well resolved by 
including the color-octet effect "J]. However, the NRQCD 
prediction[5] of the dominate transverse polarization for 
J/tp at sufficient large transverse momentum (pt) is in 
bad agreement with the almost un-polarized results mea- 
sured by CDF collaboration [1, 0]- To reveal the J/tp 
production mechanism, recently much theoretical effort 
has been made and some substantial progress has been 
achieved. At Hadron collider, the next-to-leading order 
(NLO) QCD corrections for conventional J/tp produc- 
tion gg J/'ipg in the CSM are calculated in Ref 
and the results show the pt distribution is significantly 
enhanced in large pt range; gg —^ J/tpcc in the CSM was 
found to make sizable contribution to pt distribution 9] ; 
It is reported in Ref that the polarization for the 
conventional J/tp production via the CSM turns from 
much transverse at leading-order (LO) into much longi- 
tudinal at NLO, while the J/tp produced via and 

■^S'J*' color-octet states still dominates transverse polar- 
ization at large pt at NLO[ll[; It is argued in Ref [12 
that the s-channel cut contribution may be a possible so- 
lution to the J/tp production puzzle. At e~^e~ collider, 
the NLO QCD corrections 13] to e+e~ J/tpcc and 
e+e^ J /tpgg in the CSM can well resolve the large 
discrepancies between the experimental data and LO the- 
oretical prediction. Some other attempts may be found 
in Ref.[l|- Even with these recent significant theoreti- 
cal progresses, the challenge to our understanding of the 



heavy quarkonium production mechanism is still exist. 

Since J/tp{\ ) can couple with a photon, pp{p) —> 
J* (J/tp) + X with J/tp from a virtual photon 7* frag- 
mentation can contribute to J/tp hadroproduction. It is 
thought to be small at first sight, and has never been 
taken into consideration. However, there are really a few 
examples in which the kinematic enhancement can com- 
pensate the suppression factor (a/ag)'^ that is associated 
with the QED and QCD coupling constants. For exam- 
ple, the cross section is very large in J/tp electromagnetic 
production at e+e" colliders [ij]. J/tpcc produced in 
e+e~ annihilation through two virtual photons [l^ will 
prevail over that through one virtual photon when the 
center-of-mass energy y/s > 20 GeV. And the cross sec- 
tion for e+e" — > J/tpJ/tp^^^ at LO is also a few times 
larger than that of e+e~ ^/V"7c[l2l] at v/i = 10.6 GeV, 
although its QCD corrections is negativefl^. It is easy to 
see that the pt distribution for the conventional J/tp pro- 
duction from the CSM at LO scales like 1/pt and that for 
the photon fragmentation J/tp production at LO is 1/pj, 
which may substantially enhance the contribution of the 
later at large pt region. Furthermore, the contribution of 
the QCD corrections for the conventional one gg — > J/tpg 
behaves as , so the [rric/ptf' kinematic enhancement 
of this photon fragmentation production might even com- 
pensate the suppression factor o? /ot/, compared to the 
former when pt is large enough. Therefore in this letter, 
we study the photon fragmentation J/tp production and 
polarization up to NLO in ols at the Tevatron and LHC. 
The Feynman Diagram Calculation package (FDC) pj| is 
used in the calculation. 

According to NRQCD factorization approach0|, J/tp 
production rate at hadron collider is expressed as: 

a\pp -> J/-p + ^] = X! y" dxidx2Gi/pGj/p 

•Ka[i+3^{cc)r,+X]{0^), (1) 
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FIG. 1: Typical Feynman Diagrams for J/i/> production 
through photon fragmentation at L0(B1, B2), virtual cor- 
rections(Vl,V2) and real corrections(Rl-R5). 



where p is either a proton or antiproton. The short dis- 
tance part a represents the partonic production of cc with 
quantum number n. And (C^) is the non-perturbative 
long distance matrix elements that parametrize the tran- 
sition of the cc pair into J /ip. 

At LO (a^a^), there are two partonic processes: 



Bl : 
B2 



gipi) + q{p2) 
q{vi) + q{p2) 



where q represents all possible light quark u, d, s and 
anti-quark u, d, s. The typical Feynman diagrams are 
shown in Fig.l. The partonic tree level results of the two 
processes for certain light q are trivial and given as: 
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is the 



t = 

mass of charm quark. 

The NLO correction is of a^a^ order. It contains the 
virtual(V) part and the real(R) part. At this order, the 
■^S'j^'' state could couple to one photon or one photon 
and two gluons. Only the one photon diagrams have 
kinematic enhancement at large pt region, for which the 
representative ones are displayed in Fig.l. And they also 
form a gauge- invariant subset. At present, we drop the 
contribution of the diagrams with cc pair coupling to 
three gauge bosons. The complete NLO corrections and 
their interference with full QCD processes will be dis- 
cussed in further work. 

There are 19 fragmentation diagrams among all 33 
NLO virtual correction diagrams for each Born process 
including the counter term diagrams. And the partonic 
virtual corrections for the differential cross section are: 



dt 



cx 2Re(M^Af^*) 



(3) 



The ultraviolet (UV) and infrared (IR) divergences will 
appear m the calculation of M^*. Z) = 4 - 2e dimen- 
sional regularization is applied in all the calculations. We 
adopt the same renormalization scheme as in Ref . [20| . 



The Coulomb singularity is isolated by the small rela- 
tive velocity v and then is absorbed by the long-distance 
matrix element. The IR divergences will disappear when 
including the real corrections. 

There are a few subprocesses in the real corrections. 
We divided them into seven categories below: 



99 J/tpqq, 
qq J/^qq, 

qq' J/^pqq', 
gq{q) J/ipgq{q), 



qq^ J/i^99, 
qq^ J/H'q', 
qq' J/i'qq', 



where g, q', (g, q') denote light quark(anti-quark) with dif- 
ferent flavors. For all these processes, there will be soft 
and collinear poles in the phase space integration. We 
separate them in D dimensions for coherence using the 
two-cutoff phase space slicing method [2l|. in which the 
phase space is partitioned into three parts, the soft re- 
gion, hard collinear region and hard non-collinear region 
by two small parameters: 5s the soft cutoff and 5c «C 5s 
the collinear cutoff. Then the real corrections turn to be: 



HC 



jHC 

add 



(4) 



where a^^ is the contribution of soft region including soft 
divergences and is calculated analytically in the limit of 
the soft gluon. Note that the soft poles are all associated 
with light quarks or gluons, for only the one photon frag- 
mentation diagrams are considered, a^'-^ from the hard 
collinear region contains collinear singularities which are 
also factorized out in D-dimensions. The initial state 
collinear poles are absorbed into the redefinition of the 
parton distribution function (PDF) (traditionally called 
mass factorization [2a}. Here a scale dependent PDFs 
with AIS convention [2 1| are used. After the redefinition 
of PDF, there will be an additional part cr^J^ left. The 
final state collinear poles together with the soft ones will 
cancel that of the virtual corrections, i.e. + a^'-^ + 
is finite. And the cancellation is also completed analyti- 
cally. The last part a^'~^ is finite and computed numeri- 
cally. 

To obtain the pt distribution of the Jf-ip production 
and polarization, the integral variables are transformed 
from dx2dt to Jdptdy. Then we have: 



da 
dpt 



JdxldyGa(xl,^lf)Gp{x2,^if) 



da 
dT' 



(5) 



where y is the rapidity of J/'0 in the laboratory frame 
and /i/ is the factorization scale. And the polarization 
parameter a is defined by: 



a{pt) 



daT/dpt - 2daL/dpt 
dar/dpt + 2daL/dpt 



(6) 



To calculate a{pt), the J/ip polarization vectors e(A) are 
kept through our calculation which makes the calculation 



Thereafter ct represents the corresponding partonic cross section. 
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FIG. 2: (a) The /i dependence of the cross section for the photon fragmentation J/t/; production at the Tevatron(lower) and 
LHC(upper) at LO(dashed) and NLO(sohd) with ~ iJ.f = [s.- The -pt distribution for J/?/> production at the Tevatron(b) and 
LHC(c). "QED" refers to the photon fragmentation production, "QCD" refers to the conventional production from the GSM, 
and "total" refers to the QED+QGD. The center-of-mass energies are 1.98 TeV at Tevatron and 14 TeV at LHG. 



become more complicated. And the expression of the 
polarized J/ 1/' partonic differential cross section could be 
explicitly written as: 

^ = ae(A).e*(A)+ ^ a„p, • 6(A)p, • 6*(A), (7) 

where A = Ti,T2,L. e{Ti),e{T2) and e(L) are the two 
transverse polarization vectors and the longitudinal one 
for J/^ respectively. Meanwhile, the sum over polariza- 
tions of the other particles is done in D-dimensions. One 
could find that a and j will be finite if the virtual and 
real corrections are properly dealt as mentioned before. 
To ensure the validity of our results, we also checked the 
gauge invariance by replacing the gluon polarization vec- 
tors with its momenta in the numerical calculation. 

In the numerical calculations, the CTEQ6L1 and 
CTEQ6M PDFsflil and the corresponding fitted values 
asiMz) = 0.130 and as{Mz) = 0.118 are used for LO 
and NLO predictions respectively. QED coupling con- 
stant a — 1/128 and rric = ~ 1.5GeV are chosen. 

The long distance matrix element {O'l^^) — 1.35GeV'^ is 
extracted from the leptonic decay of J/ip at QCD NLO 
level. The renormalization scale fir and factorization 
scale fif are used as /i^ = /i/ = /io = \/ (2mc)^ -|- Pt- 
The two phase space cutoffs 6s = 10^^ and Sc = Ss/50 are 
chosen, and the invariance for different values of Sg and 6c 
is obviously observed within the error tolerance. All the 
results are restricted to the NRQCD applicable domain 
Pt > 3 GeV, and \yj/^\ < 3 for the LHG, \yj/,p\ < 0.6 for 
the Tevatron. 

The numerical results for the photon fragmentation 
J/4> production are presented in Fig. 2 and 3. For com- 
parison, the plots for the conventional J/ip production 
from the GSM at NLO are also shown in the two figures, 
and they will be referred as the conventional one in the 
following statement. The curves in Fig. 2. a show that the 



scale dependence is well improved at NLO at both the 
Tevatron and LHG, and the QGD corrections enhance 
the total cross section about 30% ~ 40% for our default 
choice of scale fj,r = — Ho- From Fig. 2.b, it is clear 
that the pt distribution of J/ip production for the new 
mechanism at NLO is comparable to the conventional 
one at middle pt range {pt — 20GeV), and is larger than 
the conventional one as pt > 26GeV, and reaches about 
6 times of the conventional one when pt = 50 GeV. The 
similar results for the LHG are shown in Fig. 2.c. The 
Pt distribution for the new one at NLO is larger than 
the conventional one as pt > 35 GeV and reaches about 
10 times of the conventional one at pt = 100 GeV. The 
Pt distribution of J/ip polarization at the Tevatron and 
LHG are shown in Fig. 3. As can be seen from both 
the Tevatron and LHG results, because of J/ip from pho- 
ton fragmentation, the polarization parameter a at LO is 
positive and turn close to 1 quickly when pt increase. And 
its QGD corrections only make a slightly lower down. 
It has been given in Ref. JO] that the J/ip polarization 
at NLO is mainly longitudinal for the conventional one. 
When the contributions for the new one and conventional 
one are combined, a is negative at low pt region, and the 
photon fragmentation contribution to a becomes more 
and more important as pt increase. It makes a go from 
negative to positive gradually. 

In summary, we have suggested a new mechanism for 
J/ip production at hadron colliders, pp{p) — > ^*{J/ip) + 
X, with J/ip from photon fragmentation and calculated 
the J/ip hadroproduction from this mechanism up to 
QGD NLO at the Tevatron and LHG. Although the ob- 
tained total cross section is small, its contribution has 
large impact on the pt distribution of J/ip production 
and polarization in large pt region. The results show that 
the Pt distribution for the new production mechanism at 
NLO is larger than that of the conventional J/ip produc- 
tion from the GSM at NLO when pt > 26 (35) GeV at the 
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FIG. 3: The pt dependence of the polarization parameter a for J/tp production at Tevatron(left) and LHC (right). 



Tevatron (LHC) and reach about 6 (10) times of the con- 
ventional one when pt = 50 (100) GeV at the Tevatron 
(LHC), in spite of a suppression factor (a/cts)^. In addi- 
tion, the Pt distribution of J/ip polarization changes from 
longitudinal to transverse as pt increased by adding the 
new contribution to the conventional one. It is closer to 
the experimental result than the conventional one. In the 
future ATLAS and CMS experiments at LHC, the mea- 
sured Pt distributions of J/ip production and polarization 



can be extended to much larger pt region with its 14TeV 
center-of-mass energy and high luminosity. Therefore it 
is a important mechanism for J/ip production at large pt 
region especially for the LHC. 
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